Remarks on the shape transition from spherical to deformed gamma unstable nuclei 
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Energies and transition probabilities for low lying states in ^'''*Ba and '^"''Ru were calculated within 
a hybrid model. The ground and the first 2"*" states are described alternatively as a harmonic and 
anharmonic vibrator states while the remaining states as states with -E(5) symmetry. Results for 
^■^■'Ba are compared with those predicted by some other methods. One concludes that a gradual 
setting of the 'critical' potential yields a better agreement with the experimental data. Very good 
agreement with the data is obtained for ^''*Ru. Comparing the present results with those of E(5) 
symmetry, it is conspicuous that the present formalism add corrections to the E(5) formalism by 

- ; - bringing the predictions closer to the experimental data. Analytical relationship between the states 

^'bj^ I with (7(5) symmetry and those given by the i?(5) description is established. 

5 ■ PACS numbers: 21.10.Re, 03.65.Ge, 21.60.Fw 

CN . 

> : 

O ■ I. INTRODUCTION 

^: 

' Since the liquid drop model was developed 1], the quadrupole shape coordinates were widely used by both phe- 
nomenological and microscopic formalisms to describe the basic properties of nuclear systems. Based on these coor- 
' dinates, one defines quadrupole boson operators in terms of which model Hamiltonians and transition operators are 
^ , defined. Since the original spherical harmonic liquid drop model was able to describe only a small amount of data 
• for spherical nuclei, several improvements have been added. Thus, the Bohr-Mottelson model was generalized by 
Faessler and Grciner,2] in order to describe the small oscillations around a deformed shape which results in obtaining 
a flexible model, called vibration rotation model, suitable for the description of deformed nuclei. Later on Q this 

■ picture was extended by including anharmonicities as low order invariant polynomials in the quadrupole coordinates. 
With a suitable choice of the parameters involved in the model Hamiltonian the equipotential energy surface may 

■ exhibit several types of minima Q like spherical, deformed prolate, deformed oblate, deformed triaxial, etc. To each 
equilibrium shape, specific properties for excitation energies and electromagnetic transition probabilities show up. 
Due to this reason, one customarily says that static values of intrinsic coordinates determine a phase for the nuclear 
system. A weak point of the boson description with a complex anharmonic Hamiltonian consists of the large number 
of the structureparameters which are to be fitted. A much smaller number of parameters is used by the coherent state 
model (GSM) Q which uses a restricted collective space generated through angular momentum projection by three 
deformed orthogonal functions of coherent type. The model is able to describe in a realistic fashion transitional and 

_ ^ well deformed nuclei of various shapes including states of high and very high angular momentum. Various extensions 
^ , to include other degrees of freedom like isospin 0, single particle ,7] or octupole degrees Q of freedom have been 

■ formulated 0- 

5^ \ It has been noticed that a given nuclear phase may be associated to a certain symmetry. Hence, its properties may 
be described with the help of the irreducible representation of the respective symmetry group. Thus, the gamma 
unstable nuclei can be described by the 0(6) symmetry jlO|, the gamma triaxial nuclei by the rigid triaxial rotor 
D2 symmetry the symmetric rotor by the 5t/(3) symmetry and the spherical vibrator by the U{b) symmetry. 
Thus, even in the 50's, the symmetry properties have been grea tly appreciated. However, a big push forward was 
brought by the interacting boson approximation (IB A) |12L Il3j|. which succeeded to describe the basic properties 
of a large number of nuclei in terms of the symmetries associated to the system of quadrupole (d) and monopole 
(s) bosons which generate a C/(6) algebra. The three limiting symmetries ?7(5), 0(6), SU{'i) mentioned above, are 
dynamic symmetries for U{Q). Moreover, for each of these symmetries a specific group reduction chain provides the 
quantum numbers characterizing the states, which are suitable for a certain region of nuclei. Besides the virtue of 
unifying the group theoretical descriptions of nuclei exhibiting different symmetries, the procedure defines very simple 
reference pictures for the limiting cases. For nuclei lying close to the region characterized by a certain symmetry, the 
perturbative corrections are to be included. 

A nice classification scheme was provided by Gasten |0, who placed all nuclei on the border of a symmetry triangle. 
The vertices of this triangle symbolize the C/(5) (vibrator), 0(6) (gamma soft) and 5t/(3) (symmetric rotor), while 
the legs of the triangle denote the transitional region. Properties of nuclei lying far from vertices are difficult to be 
explained since the states have some characteristics of one vertex while some others are easy to be described by using 
the adjacent symmetry. When the mixture of the extreme symmetries is maximal, one says that the critical point 
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for the transition from one phase to another has been reached. In Ref. |l5l Il6| . it has been proved that on the 
C/(5) — 0(6) transition leg there exists a critical point for a second order phase transition while the U{5) — SU{3) leg 
has a first order phase transition. 

Recently, lachello |0, 0| pointed out that these critical points correspond to distinct symmetries, namely E{5) 
and ^(5), respectively. For the critical value of an ordering parameter, energies are given by the zeros of a Bessel 
function of half integer and irrational indices, respectively. 

It is worth mentioning that this procedure, is frequently used for quantizing the classical electromagnetic field by 
confining it inside a sphere of radius R on whose surface the solution of the field equation vanishes. Moreover, the 
magnitude of the radius determine the photon energy. 

Remarkable is the fact that the experimentalists found representatives for the two symmetries. To give an exam- 
ple, the relevant data for -^^''Ba and ^^^Sm |23| suggest that they are close to the X{5) and E{5) symmetries, 
respectively. Another candidate for E(5) symmetry, proposed by Zamfir et al. is -'^^^Pd. Using a simple IBA 
Hamiltonian, in Ref. [23| . the low lying spectrum of ^°^Pd is realistically described. Comparing the E(5) predictions 
with the experimental data concerning energy ratios in the ground band and the normalized E2 transition probabilities 
for the states 4+ and Oj, one concludes that this nucleus is a good E(5) candidate. However, in order to decide which 
Pd isotope is closer to an E(5) behavior, further investigations are necessary. A systematic search for E(5) behavior 
in nuclei has been reported in Ref.|2lj|. 

Short after the pioneering papers concerning critical poin t symmetries appeared, some other attempts have been 
performed, using other potentials like Coulomb, Kratzer |2J| and Davidson potentials These potentials yield also 
Schrodinger solvable equations and the corresponding results may be interpreted in terms of symmetry groups. 

In this paper we refer to the -E(5) symmetry and try to get a consistent description of some experimental data for 
-'^^'^Ba and ^""'Ru. The basic idea of our work consists of the fact that going from the spherical vibrator (C/(5)) to the 
0(6) like nuclei, by increasing the 'order parameter', energy levels do not reach the critical behavior at the same time 
i.e., for the same value of the above mentioned parameter. 

The fact that the critical point might be state dependent has been considered by several authors in different contexts. 
In what follows we shall mention few relevant studies, a) The dependence of the critical points on angular momentum 
and how this dependence influences the energies has been analyzed in Ref. J^J for the Davidson ,26.J potential, b) In 
Ref. [23, we analyzed the A dependence of the spectroscopic properties of Gd isotopes within the CSM approach. 
The conclusion was that the transition from spherical (SU(5)) to deformed (SU(3)) phases takes place gradually, 
i.e. not simultaneously for all energy levels. This feature has been analyzed in great details using both the CSM 
predictions and the experimental data. Therein, we stated that the critical behavior concerning the above mentioned 
phase transition is met by ^^''Gd. c) A similar conclusion concerning the dependence of the phase transition on 
the nuclear system state was reached in Ref. _28J , where some isotopes of Xe and Ba have been described within the 
CSM and the triaxial rotation- vibration model (TRVM) . Indeed, the ground state 0+ predicted by CSM behaves like 
a gamma unstable state while the state 10^ has a maximum gamma probability for 7 — 45*^. In gamma band, the 
model state for 2^^ describes a triaxial nucleus (7 — 30'') while increasing the spin, the probability distribution in the 
variable 7 goes gradually to the situation when it has two equal maxima at O'' and 60''. The 7 asymmetric features are 
described by mixing the ground and 7 model states of equal angular momenta through the anharmonicities involved 
in the boson Hamiltonian. d) The triaxial rotor Hamiltonian, with the moment of inertia depending on the nuclear 
deformation /3 in the manner predicted by the liquid drop model p9| , has eigenvalues which are functions of / (angular 
momentum) and (3. For / < 8, the energy has a minimum in /3 = 0, while for / > 8 it has a deformed minimum. 
Therefore a nuclear system described by such a Hamiltonian stays in a spherical phase if J < 8, and in a deformed 
phase for I > 8. e) A much more complex situation, when a triaxial rotator is cranked around an arbitrary axis is 
studied in Ref. |30l| . It is shown that the phase portrait depends strongly on the constraint, i.e. on angular momentum 
of the state, f) In Refs.f^ls^ it was pointed out that in nuclei exhibiting U(5) and 0(6) shape coexistence, due to 
the intruder states, the components having the two symmetries respectively, are selectively mixed, which results in 
having, for the given nucleus, states with pure symmetries and states where the two symmetries are mixed with each 
other. Very strong arguments supporting the state dependence of the critical point is also brought by microscopic 
descriptions, g) It is well known that a nuclear rotating system may exhibit a "superconducting" phase up to a 
certain critical angular momentum when a nucleon pair is broken. From this angular momentum on, one develops a 
new band of a two quasiparticle nature, having properties which differ essentially from those characterizing the states 
of lower angular momenta|3^. h) According to Ref.js^l the phase transition from a reflection symmetry to a static 
octupole deformed phase is taking place in an excited state, i) Moreover, in Ref. 35] several cases were depicted, 
where nuclear systems may stay in a phase with a good reflection symmetry for any angular momentum in the ground 
band, but at a certain angular momentum in an excited band, the static octupole deformation is set on. All examples 
mentioned above plead coherently for the idea adopted here saying that nuclear systems may change the phase when 
is promoted in an excited state. 

In order to see the effect of angular momentum dependence of the critical points on both energies and E2 transitions, 
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we suppose that the ground state ( 0+) and the first excited 2+ state are described by a spherical vibrator, while 
the remaining states correspond to a different potential energy, mocked up through a five dimensional infinite square 
well. With this assumption we calculate the excitation energies in units of the first 2+ state excitation energy as 
well as the transition probabilities of the low lying states. At the second step, we introduce some anharmonicities in 
the structure of the states 0+ and 2+. A certain mathematical relationship between the states of the U{5) and E{5) 
symmetries are pointed out. 

The aim sketched above is accomplished according to the following plan. In Section 2 the results for U(5) and E(5) 
symmetries are shortly presented. Energies and transition probabilities will be calculated in Section 3 by supposing 
that the first two states are given by a spherical vibrator while the remaining ones by the solutions of the Schroedinger 
equation for a five-dimensional infinite square well. Also, results for the situation when the ground state and the first 
excited state involve some anharmonicities, are given. In Section 4, the connection between the spherical vibrator 
states and the E(5) basis states is established. The final conclusions are drawn in Section 5. 



II. BRIEF DESCRIPTION OF THE SPHERICAL VIBRATOR AND E(5) SYMMETRIES 

Written in the intrinsic frame of reference, the original Bohr-Mottelson Hamiltonian has the expression: 



H = -- 
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/34a/3^ dl3 ' ^2 sin 37 97^"^ ^^7 4/32 ^ sin2(7 - |7rfc) 



- y - 



(1) 



where the dynamic deformation variables are denoted by (3 and 7 while the intrinsic angular momentum components 
by Qi, with i = 1,2, 3. If the potential energy term is depending only on the beta variable 



y(/3,7) = C/(/3), 



(2) 



the eigenvalue equation associated to H ^ can be separated in two parts, one equation describing the beta variable 
and the other one the gamma deformation and the Euler angles — (6*1, 62, O3). The equation in (3 is: 



A 

/32 dp^ 8(3 /32 



u{l3) 



fiP) = efiP), 



(3) 



where A is the eigenvalue of the Casimir operator of the 50(5) group. This is related with the seniority quantum 
number t, by A = t{t + 3). The 'reduced' potential u((3) and energy e are defined as: 



E 
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where E denotes the eigenvalue of the Hamiltonian H corresponding to the potential U{(3). A full description of the 
eigenstates of the Bohr-Mottelson Hamiltonian satisfying the symmetry U{5) D 5*0(5) D 50(3) D 50(2), may be 
found in Refs.jl^. In particular, the solution of the radial equation with w(/3) = P'^ is easily obtained by bringing 
first Eq. ( to the standard Schrodinger form by changing the function f to "0 by: 



V'(/3) = P^fiP). 

The equation obeyed by the new function ■(/;, is: 

d20 
d/32 

This equation is analytically solvable. The solution is: 



(r + l)(T + 2) 



/32 



e„ = 2n + T + 5/2, n ^ 0, 1, 2, r = 0, 1, 2, 3, ... 



(5) 
(6) 

(7) 

(8) 



where L'^ denotes the generalized Laguerre polynomials. The number of polynomial nodes is denoted by n and is 
related to the number of the quadrupole bosons (N) in the state, by: N = 2n + r. Consequently, the initial equation 
Q has the solution 



(9) 
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The spectrum, given by Eq. corresponds to the unitary representation of the SU{1, 1) group having the Bargman 
index fc = (r + 5/2)/2. The standard generator for SU{1, 1) are: 



(r + l)(r + 2) d , 



(10) 



where 



^2 

dp- 



^0 = 



(t + 1)(t + 2) 

/32 



The commutation relations of Hn are: 



If u{l3) ^ (3^ + gP^^, then (2.6)-(2.11) hold for r replaced [23 by 



1/2 



(11) 



(12) 



(13) 



Now, let us turn our attention to the situation considered by lachello in Ref. [T^ , where the potential term associated 
to the spherical to gamma unstable shape transition is so flat that it can be mocked up as a infinity square well 



uiP) 



0, P<Pu 

00, P > l3u 



A more convenient form for the equation in /3, is obtained through the function transformation: 
The equation for ip is 



d'^ifi 1 dip 
dp^ Pdp 



(t + 3/2) 
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p^ 



Lp = Q. 



Changing the variable /3 to z by 

z — kp, k ~ \fk 

and denoting with ^(z) — </j(/3) the function of the new variable, one arrives at: 



d'^ip 1 dip 
dz^ z dz 



lp = 0. 



(14) 



(15) 



(16) 



(17) 



(18) 



This equation is analytically solvable, the solutions being the Bessel functions of half integer order, Jt+3/2{z)- Since 
for P > Pw the function ip is equal to zero, the continuity condition requires that the solution inside the well must 
vanish for the value of p equal to pw This, in fact, yields a quantized form for the eigenvalue E. Indeed, let x^,r be 
the zeros of the Bessel function J^, : 



Jr+3/2{xi.r) =0, ^ = 1, 2, T = 0, 1, 2, ... 

Then, due to the substitution introduced in Eqs.(2.12) and (2.4) one obtains: 
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(19) 



(20) 



Concluding, the differential equation for the beta deformation corresponding to an infinite well potential provides the 
energy spectrum given by Ea. (|20|l and the wave functions: 



Pw 



(21) 



where C^.t- is a normalization factor. 
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III. A HYBRID MODEL FOR THE U(5) TO 0(6) TRANSITIONAL NUCLEI. 



Customarily, the extreme values for the energy function associated to a model Hamiltonian to be used for the 
description of the phase transition between two distinct symmetries regime is determined variationally for each state 
of the spectrum. As examples we mention the variational moment of inertia (VMI) model [3^. the Lipas-Holmberg 
equation 13 71 for the ground band energies in transitional nuclei, the solutions of liquid drop equation with Davidson 
potential j25| , cranking formalism for a many body Hamiltonian [s^ ■ In Ref . ) it is suggested that a five dimensional 
well potential, corresponding to the E(5) symmetry, may simulate a critical flat potential in the beta variable. Indeed, 
one can imagine a family of potentials, in the beta variable, depending on an order parameter rj having the following 
property. Potentials with the order parameter smaller than the critical value (rjc) have spherical minimum while 
those for rj > rjc exhibit a deformed minimum. It is remarkable the fact that the states with E(5) symmetry provide 
energy ratios and normalized B(E2) values, for a transition operator linear in /3, which do not depend on the width 
of the potential well. Note that in Ref.Iache2, the critical potential corresponding to 77c is simply mocked by the five 
dimensional well, but ignoring the process of how the order parameter is changed so to reach the value rjc- The change 
of 77 toward rjc might be caused by modifying the number of neutrons or by promoting a given nucleus in a state of 
higher angular momentum, or by changing both the neutron number and angular momentum. As a matter of fact 
the lettest version was considered in Ref. 27]. Indeed, therein we found out that for a given angular momentunm the 
critical behaviour is reached in a certain isotope while for a higher angular momentum the critical features appear in 
the next isotope. 

Restrictions required for energy extreme points to exist define in the Hamiltonian parameter space, regions charac- 
terized by the same phase space portrait for the extreme points. That means that the phase space coordinates minimal 
values have similar dependence on the structure coefficients for all points of parameter space in the selected region. 
Due to this feature the subset of the parameter space with this property is called "phase" of the nuclear system. 
Two distinct phases are separated by curves reached by changing the structure parameters. The points of separa- 
trices are called critical points. When two curves bordering the same phase intersect each other the common point 
defines the so called bifurcation point. Depending on the initial conditions the system may evolve toward one of more 
than two possible phases. A complete analysis of the nuclear phases associated to some particular phenomenological 
Hamiltonians have been presented in Refs. [4^l50ll5]| . 

It is noteworthy the fact that in order to describe several nuclear phases with a single Hamiltonian one needs to 
include anharmonicitiesfH^ . Among them the third order boson Hamiltonian plays a special role. Indeed, is this term 
the one which causes instabilities for the equation of motion for the intrinsic degrees of freedom, /3 and 7 js^ ^3 • In 
the model used by lachello |T^, there is no polynomial anharmonicity involved. By contrary, from the Bohr-Mottelson 
Hamiltonian for a harmonic liquid drop, one drops out the harmonic potential energy, which is the quadratic invariant 
/3^. As a result, one remains only with the kinetic energy operator whose eigenstates can be factorized, the radial 
factor being proportional to a Bessel function of half integer order. By choosing an infinite wall for the potential 
in beta, the boundary condition provides a quantized expression for energy. Note that such a system is quite soft 
against any deformation process indeed, since there is no surface tension term and, on the other hand, no Coulombian 
repulsion between protons. The nuclear surface is changing freely provided the deformation is not exceeding the limit 
of (3uj. The spectrum obtained in this way is quite different from the harmonic one, specific to the U(5) symmetry. 
This can be seen not only from the fact the the ratio E/^+/E2+ deviates from 2, but also from that the functions 
associated to the states 0+ and 2+ in the i?(5) and ?7(5) symmetry descriptions respectively, are quite different from 
each other. For illustration we present, in Fig. 1, the graphs of the wave functions for the two states as given by the 
corresponding £'(5) and C/(5) irreducible representations. On the other hand, differences in wave functions manifest 
themselves in the description of the transition probabilities which, as a matter of fact, define the objectives of this 
Section. 

Here we suppose that the potential energy in the beta variable is depending on angular momentum in the following 
way: 



r /32, ff < /3 < 00, L < 2, 
u{(3) = <^ 0, ff 0</3</3^, L>4, (22) 
I 00, if < 13 < 00, L> i. 



The states of interest and their energies have the following expressions: 



\LtM) = y r(n+'r + 5/2) ^^^"^'^'(^'^^"'^^'^-'(^'"^' 
^2 

Enr = — (2n + r + 5/2), (n,r) = (0,0), (0,1), L = 2t, 
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2 2 

Ei.r = (e,r) = (l,2), (1,3), (2,0). (23) 

The factor functions depending on the beta variable are solutions of Eq.© with the reduced potential given by Ea. (|22|) . 
The equation for 7 deformation and Eulerian angles (il) has the solution G^*^. Analytical solutions for this function 
for an arbitrary (nr) have been found in Ref.[3Qi] while the matrix elements for monomials in quadrupole coordinates 
between these states are given in Ref.|4l|. For low lying states, recursion formulae determining the gamma dependence 
of the wave functions were presented in Ref. ■ For the sake of completeness we listed the needed functions G^*^ in 
Appendix A. 

To calculate the B{E2) values corresponding to the electromagnetic transitions which are relevant for pointing out 
the specific features for the phase transition from U{b) to 0(6) symmetry, we need an expression for the transition 
operator. In the present work, this operator is taken as being proportional to the quadrupole collective coordinate 
which, in terms of intrinsic coordinates and Eulerian angles, has the expression: 



2M 



tf3T2Mh, n), r2M(7, ^) = (^cosjDljoin) + sin7 (Dlj^in) + i?L_2(f^))) ■ (24) 



Using the convention of Rose |43| for the Wigner Eckart theorem, the reduced transition probability between the states 
of angular momenta L and L' is obtained by squaring the corresponding reduced matrix elements of the transition 
operator: 

B{E2;L^ L')=[{L+J\T2\\l'+^)]\ 77 = 0,1,2; n = 0, 1 ; t= 1,2,3. (25) 

The states involved in the above equation are those defined in Ea. H23l) . There are some recent calculations (441] . 
which use an anharmonic structure for the transition operator. Also in the quoted paper, all states are multiplied 
by a 'condensate' factor (see Eq.(14) of the quoted reference) which contains an additional dependence on beta and 
gamma variables. The role of this condensate factor consists of transforming the boson transition operator with an 
anharmonic structure into a classical function depending nonlinearly on f3. 

The reduced matrix elements of the transition operator can be written as a product of two factors one due to 
the integration on the beta variable and one reduced matrix element of the factor depending on 7 and f2 variables, 
involved in the transition operator. The latter one can be analytically evaluated. For the cases considered here, one 
obtains:. 



{G\2\\T2\\Gl,) = {Gl,\\T2\\Gl,) = l, 



{G%\\T2\\G\^) - ^, {Gl,\\T2\\Gl,) ^ ^. (26) 

The formalism described so far was applied to ^^''Ba and ^"''Ru. We present first the predictions for ^'^'^Ba. Thus, 
results for 4 excitation energy ratios and three B{E2) ratios are presented in the first column of Table I, for /J^, = 3.1. 
How did we arrive at this value for the size of the square well potential? In order to make a suitable choice we 
analyzed first the structure of the first states corresponding to the infinite square well potential. The best way to 
do that would be to diagonalize the model Hamiltonian in the basis {\NTaJM)} with the standard notation for 
the quantum numbers: the number of quadrupole bosons (iV),the seniority (r),the missing quantum number (a), 
angular momentum (J) and its projection on the axis z (M). That can be achieved, indeed, following the procedure 
we described in Ref. (3Q] for the triaxial rotor Hamiltonian. Here we adopt another method, namely we calculate the 
overlap of the above quoted states with the lowest four vibrational states from the basis \NTaJM). From Figs. 2 and 
3 one notices that for the chosen value Pw = 3.1, the functions for the lowest 0+ and 2+ states in the two pictures, 
E{5) and U{5), are maximally overlapping, respectively. Moreover, Fig.l shows that for this value the vibrational 
states wave functions are very small although their zeros is rigorously reached only at infinity. 

Having in mind the E{5) wave functions structure revealed by Figs. 2 and 3 one may think of improving the 
agreement with the experimental data by using for the vibrational states some anharmonic components. Thus, the 
reduced energies e for the first states 0+ and 2+ are obtained by diagonalizing the anharmonic Hamiltonian: 

Hk^{N+^)+x ((6^4)0 + (6262)0) , (27) 
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-E J. 
-E + 



-E , 



S(E2;'4+,j 



B(E2;0+„-2+^) 



B(B2;2; 



B(B2;0+ -2+ ) 



B(B2;2+ 



present 



a) 



2.58 
3.71 
4.46 
0.83 

1.88 
0.28 
3.11 



present'') 



2.39 
3.47 
4.2 
0.83 

1.20 
0.55 
2.76 



*Ba (Exp.) 



2.32 
3.57 
3.66 
0.97 

1.56(18) 
0.42(12) 



E(5) 



2.20 
3.03 
3.59 
0.84 

1.68 
0.86 
2.21 



Sextic model 



2.39 
3.68 
3.70 
0.99 

1.70 
1.03 
2.12 



,4 d) 



2.09 
2.39 
3.27 
0.73 

1.82 
1.41 
2.52 



TABLE I: Results for the relevant excitation energies given in unites of the excitation energy of the first 2"*" state are listed. 
Also, the BE2 values characterizing the transitions 4^ — > 2^, Oj — » 2f and Of — » 2^ are given in units of B{E2;2f 0]^). 
"'Data from the first column are obtained by assuming that the states and 2]*" are pure vibrational states and taking for the 
states described by Bessel functions of half integer order (3^ = 3.1 In the second column the results correspond to a different 
structure for the first states. They are obtained by diagonalizing the Hamiltonian H271 with x=0.2. The other states involved 
in our calculations, correspond to (3m ~ 3.186. Experimental data, taken from Ref.,1^, are collected in the third column. In 
the fourth column are the results corresponding to the E{5) symmetry while in '^^ the fifth and sixth columns, the results of 
sextic model from Refs. |45|l and with /J"* potential |4^ respectively, are presented for comparison. 



where N denotes the number of quadrupole bosons which hereafter wiU be denoted by and 62/^- Since our 
appHcation has only an illustrative value we restrict the dimension of the diagonalization space to 2. The final results 
for energies and wave functions are: 
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+ 2x2, ^; = ^i + Ha;2^ (28) 
|0+) = Xo|0)+yo^(44)o|0), 

2V> = X2bl\0)+Y2J^{blbl)obl\0) . (29) 
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The amplitudes of the states yielded by diagonalization have the expressions 



X lu — 1 

How do the multi-boson states components, involved in the states 0+ and 2+, affect the results for energies and 
transition probabilities, can be seen in the second column of Table I. The maximal overlap of the anharmonic states 0'*' 
and 2+ with the corresponding E{5) states requires changing to the value of 3.186. Comparing the results obtained 
in this way with the corresponding experimental data one may say that anharmonicities improve the agreement with 
the experimental data. Comparing the present results with the ones obtained through different formalisms, one 
remarks that our procedure is the only one succeeding to provide a reasonable agreement with the experimental data 
for both excitation energy and the E2 decay probability of the state O2 to the first 2+ state. In this context it is 
noteworthy the fact that describing the head states of the excited bands is a decisive test for any nuclear structure 
model. 

Another nucleus whose behavior reclaims an E(5) symmetry is ^"''Ru. The even isotopes of Ru have been ex- 
perimentally studied in Ref . |47j . Their global properties which are relevant for the context of E(5) symmetry have 



Xo = 
X2 = 
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1,2 


-E + 


-B + 
1,1 

_B + 
<.3 


"to 
-E + 
"to 


E + 

"tl 

^8+ 
1,4 


-E + 
-E + 


=tl 
E + 


-E + 

"to 
-E + 
"to 


"<1 


-E I 

<0 
-E + 


E + 
^1 

°J,0 


-B + 
"to 

-E + 
"to 



'1,3 "1,0 
BCE2^4+-^2+) 



S(i32;0+'o^2+,) 



S(B2;2+j^0+^) 
S(i32;0+'3-2+,) 



B(B2;2+^^0+„) 



present 
2.40 

4.22 

6.29 

8.60 

3.49 

0.83 

1.18 
0.47 
2.85 



^Ru (Exp.) 
2.48 

4.35 

6.48 

8.69 

2.76 

0.528 

1.18(15) 
0.42(7) 



E(5) 
2.20 

3.59 

5.17 

6.93 

3.02 

0.84 

1.68 
0.86 
2.21 



TABLE 11: Results for the relevant excitation energies given in unites of the excitation energy of the first 2"*" state are listed. 
Also, the BE2 values characterizing the transitions 4^ — > 2^, Oj 2f and Of 2^ are given in units of B{E2;2'^ 0^). 
Data from the first column are obtained by assuming that the states 0^ and 2f axe eigenstates of the Hamiltonian 12711 with 
x=0.15, while the remaining ones are described by the Bessel functions of half integer order. For all states the order parameter 
has the value [3^ = 3.17 . In the second column the experimental data taken from Ref.l47l li^ are given. The results from the 
third column correspond to the E(5) symmetry. 



been analyzed within the framework of IB A model in Ref.|49|. Thus, the A dependence of the excitation energy for 
O2 suggests a transitional behavior near ^°'*Ru, although this is not reflected in the two-neutron separation energies. 
Moreover, averaging the IBA Hamiltonian on a coherent state, the resulting function of deformation and boson number 
admits N=8, which corresponds to ^"^^Ru, as a critical value. Indeed, the energy function has a minimum in /3 = for 
N < 8, while for > 8 the minimum value is reached for a nonvanishing f]. The IBA description uses a 5 parameter 
boson Hamiltonian and an anharmonic structure for the transition operator. Four of the five parameters are fixed so 
that an overall agreement for all ^4Ru isotopes is obtained. Predictions for energies and transition probabilities are 
in very good agreement with the corresponding experimental data. 

Results of our calculations for ^"^"^Ru are collected in Table II. Here we considered for the first two states, 0+ and 
2+, an anharmonic structure by diagonalizing the boson Hamiltonian (|27|l . The higher excited states are described 
by the E(5) formalism. As in the previous case there are two parameters involved, Pw and x. The first one is fixed 
so that the overlap of vibrational states and E(5) states have a maximal overlap. This way of fixing /3m is justified 
by the fact that at critical point the specific features of the extreme symmetries, SU{5) and 5*0(6), are maximally 
mixed. It is worth noting that by contrast to the E(5) formalism where the normalized energies do not depend on (3^, 
here the results for both transition probabilities and energies are depending on this parameter. The other parameter 
(x) was fixed so that the most sensitive transition probability is reproduced. For comparison in Table II we give also 
the results obtained with the E(5) formalism. Note that excepting the energy of all other data are fairly well 
reproduced. Moreover, from Table II it is conspicuous that the present results are much closer to the experimental 
data than those provided by the E(5) description. 

Concluding, the features mentioned above support our hypothesis asserting that the critical value for the Hamilto- 
nian order parameters (which is is angular momentum dependent. 



IV. RELATIONSHIP OF SPHERICAL VIBRATOR AND E(5) STATES. 

As we stated already, the infinite 5 dimensional well is simulating the situation met for some transitional nuclei of 
a fiat potential in the beta variable. This is, indeed, a very particular picture taking into account that the 0(6) nuclei 
are soft only in the 7 variable. If that is the case, then one suspects that there are some mathematical features which 
vindicate this approximation. 

Also we note that the procedure used for obtaining the quantized energy for U (5) and E{5) symmetries are somehow 



9 



different. While for U{5), one requires that the wave function has a regular behavior at infinity, for the E{5) states 
the condition is to vanish for the wall coordinate. 

Energies provided by the zeroes of the Bessel functions of half integer order are organized in bands labeled by 
the quantum number ^, which numbers the above cited zeros. For example, the ground band comprises states with 
L = 2r, and = 1. On the other hand the spherical vibrator states are classified in bands according to a different rule, 
namely the ground band is the set of n = states of highest seniority, the beta band is the second highest seniority 
band, etc. According to these schemes of classification one may put in correspondence the states of spherical vibrator 
ground band with those of the ground band in the E(5) picture. Thus, to the highest seniority ^ = 1 states would 
correspond the highest seniority n = states. Such a correspondence is also supported by Fig. 1 showing that the 
corresponding states have similar allure, although they differ in magnitude. 

Concerning the ground band energies for the spherical vibrator another remark is worth to be mentioned. Since 
the polynomial L^^'^^^(.t) is linear in x, it has only one zero at 

which can describe the ground band energies by considering 

EOt — 00- 

with /?o denoting the zero given by Eq.(4.1) where x is replaced by /3q. Inspecting Eq.Q, one sees that for these 
wave function zeros, the second order derivative is also vanishing. Reversely, the solution of the equation (|32|l in /3o 
is a zero for the wave function. This reflects a classical property of this potential consisting in that the solutions for 
Ea. (|32|l provides the turning points of the classical trajectory. The comments from above raise the question whether 
it is possible that all energy levels from the spherical vibrator picture can be given as zeros of some generalized 
Laguerre polynomials of a contracted variable. Moreover, is it possible to describe the vibrator states in terms of 
Bessel functions of half integer order? The answers of these questions are given by using the asymptotic behavior of 
generalized Laguerre polynomials |53j which approximate quite well the Bessel functions according to the formula. 



(31) 



(32) 



L{s, n, z) 
L{s, n, z) 



r(n + s 
Js{z). 



exp 



4 + 8n + 4s 



2 + 4n + 2s 



(33) 



The function L approximates exceptionally well the Bessel function of half integer for large values of n. The 
approximation is, however, reasonable good even for small values of n. For Illustration we present in Table III, the 
values of the deviations from unity of the ratio of L{s, n, z) and Js{z) for few values of n. 



R{s,n,z)~ L{s,n, z) / Js{z) — 1 



(34) 



Actually this equation fixes the conditions under which the E(5) energies and functions might be described in terms 
of U(5) energies and states, respectively. 



n 


r + 3/2 




3/2 


5/2 


7/2 


9/2 


11/2 


13/2 


15/2 


17/2 


19/2 


21/2 


23/2 





26.4 


9.8 


4.7 


2.6 


1.6 


1.1 


0.7 


0.5 


0.4 


0.3 


0.2 


1 


8.1 


3.9 


2.2 


1.4 


0.9 


0.7 


0.5 


0.4 


0.3 


0.2 


0.2 


2 


3.9 


2.1 


1.3 


0.9 


0.6 


0.5 


0.3 


0.3 


0.2 


0.2 


0.1 


3 


2.2 


1.3 


0.9 


0.6 


0.4 


0.3 


0.3 


0.2 


0.2 


0.1 


0.1 


4 


1.5 


0.9 


0.6 


0.4 


0.3 


0.2 


0.2 


0.2 


0.1 


0.1 


0.09 


5 


1.0 


0.7 


0.5 


0.3 


0.3 


0.2 


0.2 


0.1 


0.1 


0.09 


0.07 



TABLE III: The function R{s, n, z) given by Eg. 1341 . multiplied by a factor of 10^, was calculated for s = r + 3/2, z = 2 and 
several values of the quantum number n . 

Let us replace z — ^fkfi in Eq. H33|l and denote by ye^,n,r the ^-th zero of the generalized Laguerre polynomial 
^T+3/2j-^-j^ Then the corresponding energy is: 



2B 



{An + 2r + 5) 



(35) 
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and the normalized excitation energies 



Rn 



El, 



El. 



El, 



El, 



(36) 



are readily calculated. For large values of n this ratio must converge to the ratio determined by the energies given by 
the zeros of the Bessel function of half integer order. 



Rr 



no 



(37) 



no 



In Table IV we list the values of ratio (|36() and 1)3 7() . From there, one may see that although the convergence is reached 
for large values of n, even the low values of n produce energy ratios which are close to the E(5) limit. 



>^T+3/2 
jT+3/2 

7-T+3/2 

rT+3/2 
^3 

LI 



+3/2 
+3/2 
rT+3/2 



-'^lOO 



+3/2 



2.19859 

2.25 
2.23793 
2.22819 
2.22136 
2.21657 
2.21310 
2.19872 



3.58982 

3.75 
3.71562 
3.68612 
3.66488 
3.64966 
3.63852 
5.59027 



5.16941 

5.5 
5.43480 
5.37567 
5.33204 
5.30023 
5.27666 
5.17046 



6.93412 

7.5 
7.39708 
7.29876 
7.22458 
7.16969 

7.128 
6.93615 
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8.8814 

9.75 
9.60389 
9.45723 
9.34431 
9.25963 
9.19550 
8.88490 



TABLE IV: The ratios Rn.T given by Eg.ira are calculated for several values on (n,r). In the first column are given the 
generalized Laguerre polynomials whose first zero determine the involved energies. In the first row the i?(5) limit, i.e. the ratio 
Rt given by Eq. 11371 , is listed. The corresponding energies are zeros for the Bessel function given in the first column. 



Concluding, although the generalized Laguerre polynomials and Bessel functions of half integer order account for 
different symmetries, according to Eq. (|33|l they approximate each other for large values of n. However, in virtue 
of Eq.(2.11) the generalized Laguerre polynomials are eigenstate of one of the SU{1,1) generators. On the other 
hand this group may be contracted to an Euclidean one. Therefore, the asymptotic behavior reflected in Ea. H33l) is 
supported by the group relationship mentioned above. On the other hand the fact that asymptotically the generalized 
Laguerre polynomials can be approximated by Bessel functions of half integer order is a nice justification of replacing 
the potential in (3 characterizing the transitional U(5)-0(6) nuclei, by a five dimensional infinite well potential. 



V. SUMMARY AND CONCLUSIONS 



The main results obtained in the present paper can be summarized as follows. We suppose that in the critical point 
of the C/(5) — 0(6) transition the potential in (3 is angular momentum dependent and given by Ea. (|22|l . To be more 
concrete, the ground and the first 2+ states are described by a spherical vibrator Hamiltonian while the higher states 
by a Hamiltonian involving a five dimensional infinite well potential in f3. This simulates a scenario where the higher 
states undergo the phase transition while the first states still behave as spherical vibrator states. Thus, the assigned 
wave functions are those from Eq. H23|l . At the next stage we considered some anharmonicities in the structure of the 
spherical states. The results obtained for energies and transition probabilities are collected in Table I for ^'^^Ba and 
in Table II for io^r^. 

From Table I one notices a good description for the state O^q. This supports our hypothesis for a gradual transition 
from a spherical to a i?(5) symmetry. The predictions shown in Table II for ^""^Ru suggest a very good agreement 
with the experimental data. Moreover, results of our calculations are much closer to the experimental data than those 
obtained with the E(5) symmetry formalism. The present results encourage us for searching for new nuclei whose 
behavior is compatible with an E(5) symmetry in some states but with U(5) and possible with 0(6) symmetries in 
the remaining states. 

The connection between the basis states specific to the two symmetries was discussed in Section 4. The relation 
(|33|l . which holds in the asymptotic region for the quantum number n may be considered as a mathematical frame 
for the phase transition and moreover suggests a possible way of relating the descriptions of the two phases. 
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VI. APPENDIX A 



Here we list the normalized functions G„t- which we need in Section 3. 



1 



1. /5 
47rV 2 

J_ /5 
47rV 2 
1 

8ttV2 



47r' 



cos 87, 



47r V 2 
sin 7 ^^2 /ON , n2 



cos 272^2^0 (r!)- 



V2 

sin 27 ^ 2 



V2 



15 



(6cos2 7 + sin2 7) i?4„(r!) + sin 27 (Z^ial^^) + ^M,-2(f^)) 



-— Sin2 7(i?i,4(f7)+i?i,_„4(f7)) 
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FIG. 1: (Color online) First row:the wave functions describing the first two states, 0^-left panel and 2^-riglit panel, in a 5 
dimensional infinite well. Second rowithe wave functions describing the first two states, 0^-left panel and 2^-right panel, in a 
5 dimensional oscillator well. 
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FIG. 2: (Color online)The overlaps of the functions /j,^ with ^ = l,r = given by Eq.(2.20) and F„r (Eq. (2.9)) with 
(n, r)=(0,0) (upper-left panel), (l,0)(upper-right panel), (2,0) (bottom-left panel),(3,0)(bottom-right panel). 




FIG. 3: (Color online) The overlaps of the functions /j.t with ^ = l,r = 1 given by Eq. (2.20) and FnT (2.9) with 
(n, r)=(0,l) (upper-left panel), (l,l)(upper-right panel),(2,l)(bottom-left panel),(3,l)(bottom-right panel) 



